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Abstract 

We use the theory of dynamical invariants to yield a simple derivation of noncyclic 
analogues of the Abelian and non-Abelian geometric phases. This derivation relies 
only on the principle of gauge invariance and elucidates the existing definitions of the 
Abelian noncyclic geometric phase. We also discuss the adiabatic limit of the noncyclic 
geometric phase and compute the adiabatic non-Abelian noncyclic geometric phase 
for a spin 1 magnetic (or electric) quadrupole interacting with a precessing magnetic 
(electric) field. 

PACS numbers: 03.65.Bz 



1 Introduction 

Since the publication of Berry's seminal paper [0 on the adiabatic geometric phase, the 
concept of geometric phase has been generalized in a number of ways. Following the work 
of Aharonov and Anandan on nonadiabatic geometric phase, Samuel and Bhandari || 
showed that one can indeed define an analogue of the Abelian geometric phase for a quantum 
state that does not undergo a cyclic evolution. Zak ||, Aitchison and Wanelik ||, Mukunda 
and Simon |J, Pati J7L El, and de Polavieja and Sjoqvist |J have ellaborated on the theoretical 



aspects of noncyclic geometric phases, and Wu and Li [10 , Weinfurter and Badurek [ I I 



Christian and Shimony [12|, Wagh and Rakhecha [13| and Wagh, Rakhecha, Fischer, and 
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Ioffe |14| have explored its experimental consequences. In all these investigations the authors 
consider Abelian noncyclic geometric phases. The main purpose of the present paper is to 
offer an alternative approach to noncyclic geometric phases which clarifies the existing results 
on the Abelian noncyclic geometric phase and allows for its non-Abelian generalization. 

The cyclic geometric phase can be conveniently discussed within the framework of the 
theory of dynamical invariants of Lewis and Riesenfeld |T5|] . The application of dynamical 
invariants in the study of the cyclic geometric phases has been considered by Morales JTE| 



and Monteoliva, Korsch and Nunes |T7| for the Abelian case and by the present author |18 
for the non-Abelian case. 

In this article, we shall first present a brief review of the necessary results from the theory 
of dynamical invariants and comment on their relevance to the cyclic geometric phases in 
section 2. In section 3, we outline an alternative approach to nonadiabatic cyclic geometric 
phase which is essentially the same as Berry's approach to the adiabatic geometric phase 



]| . This generalizes the analysis of Ref . [19| in which such an approach is developed for the 
description of the nonadiabatic cyclic geometric phase of a magnetic dipole interacting with a 
precessing magnetic field. In section 4, we derive an expression for the evolution operator and 
discuss its gauge invariance. In section 5, we give our definition of the noncyclic geometric 
phase and explore its relationship to the cyclic geometric phase [F21], ITH]. In section 6, we 
restrict to the Abelian case and compare our definition of the noncyclic Abelian geometric 
phase with the earlier definitions In section 7, we present a discussion of the adiabatic 
approximation and show that our analysis reproduces the results of Ref. |§ in the adiabatic 
limit. In section 8, we calculate the adiabatic non-Abelian noncyclic geometric phase for 
a spin 1 magnetic (or electric) quadrupole interacting with a precessing magnetic (resp. 
electric) field. Finally, we present our concluding remarks in section 9. 

2 Invariant Operators and Cyclic Geometric Phase 



As shown by Lewis and Riesenfeld |L5| , one can write the solution of the Schrodinger equation 



ij t \m) = mm)) 



for a time-dependent Hamiltonian H(t), as a linear combination of certain eigenvectors of a 
Hermitian dynamical invariant. The latter is a Hermitian operator I(t) satisfying 

dl(t) 



dt 



i[I(t),H(t)} . (2) 



We shall assume that both H(t) and I(t) have discrete spectra. 

Now let us label the eigenvalues of I(t) by A n and the degree of degeneracy of A n by l n . 
Furthermore, let |A n , a;t) be arbitrary orthonormal eigenvectors of I(t) satisfying 

I(t)\X n , a;t) = A„|A n , a;t) , (3) 
(A m , b; t\X n , a; t) = S mn 5ba , (4) 

In 

J2J2\X n ,a;t)(X n ,a;t\ = 1 , (5) 

n a=l 

where a and b are degeneracy labels taking their values in {1, 2, • • • , l n }. 

Clearly unlike the eigenvalues A n and the corresponding degeneracy subspaces 7i\ n {t), 
the eigenvectors |A n , a; t) are not uniquely determined by the eigenvalue equation @. They 
are only determined up to unitary transformations of the degeneracy subspaces T~L\ n (t), 

In 

|A n , a; t) \X n , a; t)' = J2\X n , b; t)u ba (t) , (6) 

b=l 

where u a b(t) are the entries of an arbitrary unitary l n x l n matrix u(t). 

The main result of Lewis and Riesenfeld |T5j is that one can choose a particular set of 



eigenvectors |A n , a; t)' that are solutions of the Schrodinger equation ([!]). These eigenvectors 
are given by 

In 

\X n ,a;t)' := ^ |A n , a; t)< fe (t) , where (7) 



6=1 

u n {t) := Te- i /o An (*') d '' M "(0) = Te- i /o[ £,l ( i ')-' 4n ( i ')] d '' M -(0) 



T stands for the time-ordering operator and A n (t), £ n (t) and A n (t) are Hermitian l n x l n 
matrices with entries 

K,{t) ■= £m-A n ab (t), (9) 



^b(t) ■= (X n ,a;t\H(t)\X n ,b;t) , and (10) 

d 

di^ 



■Kbit) := i(X n ,a;t\^-\X n ,b;t) , (11 



respectively. In other words, u n (t) is a solution of the matrix Schrodinger equation 

ij t u n {t) = A n (t)u n (t) . (12) 
If the matrices S n (t) and A n (t) commute, then we can write 

u n (t) = Te-^o £n{t ' )dt 'Te- l I> n ^ dt 'u n (0) . (13) 



For the case where the invariant I(t) is periodic Jl7|, |T^], i.e., I(T) = 1(0) for some T 



\X n ,a;T) = \X n , a; 0) and \X n ,a;t)' undergoes a cyclic evolution. The corresponding non- 
Abelian nonadiabatic cyclic geometric phase [3(J is given by T n (T) where T n (t) is defined to 
be the unique solution of the matrix Schrodinger equation 

ijT n (t) = -A n (t)T n (t) , r n (0) = 1 . (14) 

Alternatively, 

T n (t) :=Te i ti An{t ' )dt ' . (15) 
If the eigenvalue A ra is nondegenerate, then l n = 1 and we have 

u n (t) = e i5n{t) T n (t) , where (16) 

S n {t) := - fs n (f)dt' = - f\\ n -t'\H(t')\\ n -t')dt' , (17) 
Jo Jo 

T n (t) = e^W, and ln {t) := /* A n (t')dt' = f i(X n] t'\^-\\ n ; t')dt' . (18) 

Jo Jo at' 

In this case T n (T) is the Abelian nonadiabatic cyclic geometric phase ||. 

3 An Alternative Approach to Nonadiabatic Cyclic Ge- 
ometric Phase 

In order to make the geometric character of the cyclic geometric phase more transparent, 
we shall express the invariant I(t) as a linear combination of a set of linearly independent 
constant Hermitian operators Xi, 

N 

!(*)=£**(*)**■ (19) 

i=l 

Here iV is a fixed nonnegative integer, the coefficients 9 l are real-valued functions, and Xi 
are generators of the group U(TC) of unitary transformations of the Hilbert space 7i. If the 



system has a finite-dimensional dynamical group G, then X( are the representations of the 
generators of G. In this case N is just the dimension of G. However, if the Hilbert space 
is infinite-dimensional, then in principle one will need an infinite number N of generators 
Xi of U{7i) to satisfy ([l9|) p], and one must find a way to make the right-hand side of (19) 



well-defined. We shall not be concerned with subtleties of the infinite-dimensional unitary 
group PHI , and assume that N is finite. 

Now since the time-dependence of I(t) is governed by those of the parameters 9 l , we 
can consider the parameter-dependent operator I [9] with eigenvalues^ A n and eigenvectors 
\X n , a; 9), and write 

I(t)=I[9(t)} and |A n , a; t) = \X n , a; 9(t)) . (20) 

Here 9 stands for (9 l , 9 2 , ■ ■ ■ , 9 N ). The parameters 9 l may be viewed as local coordinates of 
a parameter space Ai- As time progresses they trace a curve C in Ai. 

If the system possesses a dynamical group G, then we can introduce the parameter- 
dependent Hamiltonian 

TV 

H[R] = Y,R i Xi, (21) 
i=i 

and identify H(t) with H[R(t)], [|TJ. This means that the parameter space Ai of the invariant 
(|19|) is the same as the parameter space of the Hamiltonian ( 21f) . In this case, we can use 



the results of Ref. |22[ to identify Ai with a submanifold of the flag manifold G/T where T 
is a maximal torus of G. 

Now suppose that the curve C traced by the parameters 9(t) lies in a local coordinate 
patch of the parameter space Ai. In this case, we can introduce the nonadiabatic analogue 



of the non-Abelian Berry connection one-form f23| , |19 

N 



A n [9] = YsKWW ■ ( 22 ) 



i=i 



The matrix elements of A n [9] and its components A™ [9] are given by 



A n ab [9] := i{\ n ,a;9\d\\ n ,b;9) , and (23) 
(Am) ab := i(\ n ,a;9\-^\\ n ,b;9) , (24) 



1 Note that for each value of t there is a finite number TV of generators Xi for which (|l^) holds. However, 
for an infinite-dimensional Hilbert space, as t changes N might not have an upper bound. Therefore, in order 
to satisfy ( jTsj ) one would in general need to include an infinite number of generators of the group UiTL). 

2 Note that the eigenvalues of an invariant operator are constant [|l5|, Q. 



respectively. In Eq. ( p3|) d stands for the exterior derivative with respect to Q % . In view of 
Eqs. (P), (m, pa), (m, and we have 



iV rlfttt) 

A n {t)dt = Y J ^W)[-j^-dt = A n [6{t)}, and (25) 
i=i c ^ 

r n (t) = Ve^m Anm , (26) 

where P stands for the path-ordering operator. In particular, for a periodic invariant, the 
curve C traced by 6{t) is closed, and the non-Abelian cyclic geometric phase takes the form 

r n (T) =Ve i $c An W . (27) 



As pointed out in Ref. [|T8 ] , this expression agrees with Anandan's definition of non-Abelian 



cyclic geometric phase |2C 



If C does not lie in a single coordinate patch of M, one must evaluate the path-ordered 
integrals in (|26|) and (p7|) along the segments of C belonging to different coordinate patches 
and multiply the resulting unitary matrices in the order in which the curve C is traversed in 
time. 

4 Evolution Operator and Its Gauge Invariance 

In general, we can write any solution of the Schrodinger equation ([J) as a linear combination 
of \X n ,b;t)', i.e., 

W)> = EE*.M>' = E E cy n ba (t)\x n ,b;t) , (28) 

n a=l " a,b=l 

where 

C'a n :=E<a t (0)(An,&;#(0)), and C n a := (A n , a; #(0)) . (29) 

6=1 

In view of Eq. (p8|), the evolution operator is given by 



=E E wS,(t)|A n ,a;t)(A n ,6;0| . (30) 

i a,b=l 

Now let us recall that the eigenvectors \X n ,a;8) are not uniquely determined by the 
eigenvalue equation 

I[9]\\ l ,a;9) = X n \X n ,a;9) . (31) 



They are subject to arbitrary gauge transformations 

|A n , a; 9) -> |A n , a; 9)~ := £ \ X n, &; t)v&[9] , (32) 

6=1 

where v™ b [9] are entries of an l n x l n unitary matrix v n [9]. In fact, if we denote the local 
coordinate patch corresponding to the coordinates 9 % by O, v n may be viewed as a smooth 
function mapping O to the unitary group U (l n ) . 

Using Eq. (|23|), we can easily derive the gauge transformation law for A n [0], namely 

A n [9] -> A n [9] = v n [9} ] A n [9}v n [9] + iv n [6]Uv n [6] . (33) 



Furthermore, in view of Eqs. (EST), (H), and ( P5[) , we have the following transformation rules 
for T n (t) and u n {t). 

T n (t) -> f"(t) = w n ^(t)] t r n (t)w"[^(0)] , (34) 
w"(t) -> w n (t) = v n [9(t)]^u n (t)v n [9(0)} . (35) 

A simple consequence of Eqs. fl32[) and (|35|) is that the evolution operator (|30D is invariant 
under the gauge transformations. 

For a periodic invariant I(t), with 9{T) = 9(0), T n (T) is related to T n (T) by a similarity 
transformation 

f"(T) = w n [0(O)] t r n (T)t;"^(O)] . (36) 

In other words, under a gauge transformation (|32"D r™(T) transforms covariantly. Therefore, 
its eigenvalues and in particular its trace are gauge- invariant. These are essentially the 
physically observable quantities associated with the non-Abelian cyclic geometric phase. 



5 Noncyclic Geometric Phase 

The main reason for Berry's consideration of a periodic Hamiltonian Q and Aharonov and 
Anandan's consideration of cyclic evolutions 0] is the fact that for a cyclic state with period 
T- in our approach a T-periodic dynamical invariant - the unitary matrix T n (T) transforms 
convariantly under a gauge transformation fl32|) . This property of T n (T) guarantees that 
its eigenvalues and its trace are gauge- invariant. This in turn raises the issue of exploring 



their physical consequences. There is, however, another way of constructing gauge-invariant 
quantities using T n (t). 

In view of Eqs. (|34]) and (|35|) , T n (t) and u n (t) have the same gauge transformation 
properties. This means that if we replace u n (t) in the expression (|30|) for the evolution 
operator by T n (t), then we shall still obtain a gauge-invariant operator, namely 

V(t) :=EE T n ab (t)\\ n ,a;6(t))(\ n ,b;6(0)\ . (37) 

n a, 6=1 

In fact, because the gauge transformations (|3^) are the unitary transformations of the de- 
generacy subspaces 7~L\ n [0{t)] ) the restriction (or projection) of V{t) onto the degeneracy 
subspaces, i.e., 

V n (t) := £ r n ab (t)\X n ,a;e(t))(X n ,b;e(0)\ , (38) 

a,b=l 

will also be gauge- invariant. By construction, the operators V n (t) are uniquely determined 
by the curve C and its end points. In particular, they are independent of the duration of the 
evolution. Therefore, they are also geometric quantities. 

Next let us note that the gauge-invariance and geometric character of V n (t) will not be 
affected, if we exchange the positions of |A n , a; 9(t)) and (A„, b; 6*(0)| in Eq. (|3"8]). In this way 
we obtain a set of gauge invariant scalars: 

rP(t) := £ r: b (t)(\ n ,b;9(0)\\ n ,a;9(t)) = £ <(t)C(t) = trace [w n (t)T n (*)] , (39) 

a,b=l a, 6=1 

where we have introduced the Z„ x Z n matrices w n {t) with entries 

<(*):=(A n ,6;e(0)|A n ,a;^)). (40) 

By definition w n (t) only depend on the end points of the curve C. T n (t) are also uniquely 
determined by C. Therefore, as expected U n (t) are geometric quantities. 

Now let us consider a T-periodic invariant (a cyclic evolution), for which |A„, a; 0(T)) = 
|A„, a; #(0)). In this case, w n (T) = w n (0) is just the identity matrix, and Eq. (|39| ) reduces to 

IT(T) = trace[r n (T)] . (41) 

Therefore, for a cyclic evolution U n (T) yields the trace of the non-Abelian cyclic geometric 



phase p3, |4 



On the other hand, since U n (t) is gauge-invariant, we can compute it in a basis {| A n , a; 0(t))+} 
in which T n (t) is diagonal. If we denote the eigenvalues of T n (t) by e* 7 ™^, then we have 

n n (t) ,(A n ,a;0(t)|A n ,a;0(O)>* . (42) 

0=1 



As seen from Eq. (142]) , in general IT n (t) is a nonunimodular complex number. 



In view of the above analysis, the matrix f n (C) defined by 

f n (C) = t n (t) := w n (t)T n (t) , (43) 

is a gauge-covariant geometric quantity. We shall therefore identify it as the non-Abelian 
noncyclic geometric phase factor. Clearly, the eigenvalues of T n (C) and its trace namely IP(t) 
are gauge-invariant quantities which can, in principle, be observed experimentally. The non- 
Abelian noncyclic geometric phase factor T"(C) is therefore as physically significant as its 
cyclic counterpart (|27|) . 

As mentioned above, for a cyclic evolution where the invariant I(t) is T-periodic, |A n , a; 6{T)) 
| A n , a; 9(0)) and w n {T) is the identity operator. In this case, Y n {T) is identical with the cyclic 
geometric phase factor T n (T) given by Eq. (f2~7|). 

6 Abelian Noncyclic Geometric Phase 

For a nondegenerate eigenvalue A n of the invariant I(t), we have 

f n (t) = w n (t)T n (t) = U n (t) = (A n ; 0(O)|An; 0(t))e**>® , (44) 

where 7 n (t) is the phase angle given by (|i~8l) . In particular, 

|f"(t)i = K(t)i = |(A n ;0(O)|An;^))| (45) 

depends only on the end points of the curve C. If |T n (t)| 7^ 0, then we can consider the phase 
of T n (t) which is given by 

f n (t) 



e i7n(t) ;= = i[u»(t)+7»(t)] where ( 46 ) 

|r n (t)| ' v 7 

iVn{ t) ™ n (t) (x n ,e(o)\x n ,e(t)) \(\ n ,6(o)\\ n ,6(t)y 12 

^ |t^(t)| \(\ n ,6(0)\\ n ,6(t)}\ 



(47) 



(A n ;^(t)|A„^(0)) 

The phase angle j n (t) is a real noncyclic geometric phase angle. It consists of two pieces: 
7„(t) that depends on the curve C, and r] n (t) that depends on the end points of C. 



The phase factor (|46|) coincides with the 'real noncyclic geometric phase' introduced by 
Aitchison and Wanelik ||. As discussed by Aitchison and Wanelik it is equivalent to the 
noncyclic geometric phase of Samuel and Bhandari || and Mukunda and Simon || . 

7 Adiabatic Approximation and the Noncyclic Geo- 
metric Phase in the Adiabatic Limit 

Let H [R] be a parameter-dependent Hamiltonian with a discrete spectrum. We shall denote 
the eigenvalues of H[R] by E n [R], their degree of degeneracy by TV, and the corresponding 
degeneracy subspaces by TC n [R]. Let \n, a; R) form a complete set of orthonormal eigenvectors 
of H[R\. They satisfy, 

H[R)\n,a;R) = E n [R)\n, a; R) , (48) 

where a is a degeneracy label taking its value in {1, 2, ■ • • ,7V}. 

Now consider the time-dependent Hamiltonian H(t) := H[R(t)], where the parameters 
R(t) = (-R 1 (t), R 2 (t), ■ ■ ■ , R d (t)) trace a curve C in the parameter space M of the Hamilto- 
nian. We shall denote the duration of the evolution of the system by r and suppose that 
during the evolution no level- crossings occur. Furthermore, let I(t) be a dynamical invari- 
ant satisfying (H), and suppose that it depends on a set of parameters 6 = (9 1 , 9 2 , ■ ■ ■ , 9 N ), 
i.e., I(t) = l[0(t)} where 8(t) traces a curve C in the parameter space M. of the invariant. 
Since lit) yields the solution of the Schrodinger equation, the dynamics of the system can 
be encoded in the definition of a function 

F:M^M } defined by F(R):=0. (49) 

In particular, F maps the curve C onto the curve C, and I(t) = l[0(t)] = I[F(R(t))]. Note, 
however, that there does not exist a universal function F describing all possible dynamical 
processes. In other words, the definition of F also depends on the choice of the Hamiltonian 
or alternatively the curve C. In this sense, it is more appropriate to define the function 

T-Vm^Vm by T[C) = C. (50) 

where Vm and Vm are the space of paths in M and M. respectively. The latter are infinite- 
dimensional spaces. Therefore, it is more convenient to use the function F with the provision 
of its nonuniversal character. 



Next we define a normalized time variable by s := t/r, and assume that for sufficiently 
large values of r we can expand I(t) and H(t) in powers of r _1 , i.e., 

oo oo 

I(t) = I(ts) = I (s) + r' £ Ii(s) and H(t) = H(rs) = H (s) + ]T r^H^s) , (51) 

where If and Hi are Hermitian operators and Iq 7^ 7^ H 0. If we substitute t = rs and (|5l| ) 
in (0) and take the limit r — > 00, we find [I (s) , H (s)] = 0. This means that in this limit, 
where I(t) — > Iq(s) and H{t) — > H (s), the eigenvectors of /(t) and -ff(t) coincide. Using 
the notation of the preceding sections, we have 

lim \A n ,a:t) = \n,a;t) . (52) 

T— >0O 

Because |A n ,a;t) = \X n ,a;9(t)), \n,a;t) = \n,a; R(t)), and (|52D is independent of the form 
of the curve C, 

lim \X n ,a;9)\ c = \n,a;R)\ c . (53) 

In particular, this implies that for sufficiently large values of r, we can choose an invariant 
whose parameter space is the same as that of the Hamiltonian, M. = M. Therefore, C and 
C belong to the same parameter space M. In this case, we can also express (^) and (|53|) by 

lim F = id M • (54) 



where idjvf is the identity function on M. 

Now let us use Eqs. (|30D and (^) to compute the evolution operator in the limit r — ► 00. 
This leads to 

lim U(t) = U (0) (t) , where (55) 

T^OO 

TV 

U (0) (t) :=E E u: ab (t)\n,a;t)(n,b;0\, (56) 

n a,b=l 

u »(t) = e*" w r^(t) , (57) 

Son(t) ■= - f E n {t')dt> , (58) 
Jo 

TZ(t) :=Ve^S A " lR] , and (59) 

(A n [R]) ab :=i(n,a-R\d\n,b-R) . (60) 

Here the subscript o is inserted to mean that the corresponding quantities are obtained in 
the adiabatic limit (r — > 00) from the ones without a subscript o. The matrix of one-forms 
is the non-Abelian Berry connection one- form [^3[ . 



In practice, r is a finite number and the limit r — > oo is interpreted by the condition 
that t must be much larger than the time (inverse of energy) scale of the quantum system. 
If this happens to be the case the above results may be used. It is not difficult to check 
that an operator I(t) with the same eigenvectors as the Hamiltonian satisfies Eq. (|2|) only if 
the eigenvectors of the Hamiltonian are constant. This is often not the case. Therefore, in 
an eigenbasis {\n,a;t)} of the Hamiltonian, the invariant I(t) is not diagonal. However, if 
the above adiabaticity condition is fulfilled, i.e., r is much larger than the time scale of the 
problem, then the off-diagonal matrix elements of I(t) are much smaller than its diagonal 
matrix elements. The approximation scheme in which one neglects the off-diagonal matrix 
elements of I(t) is called the adiabatic approximation, f2"5] , pUj , |57|| . In this approximation, 
we have 

\X n ,a;t) « \n,a;t) , U(t) « U {0) (t) , and F « id M . (61) 
This is a valid approximation scheme if and only if 

AZ n (t) :=i(m,b;R\±\n,a;R) = * for m^n. (62) 

Here m and n are arbitrary labels (satsifying m ^ n) and a and b are arbitrary degeneracy 
labels associated with the eigenvalues E n {t) and E m (t), respectively. The second equation 



in (|62|) is obtained by differentiating both sides of Eq. (f48|) with respect to time and taking 
the inner product of both sides of the resulting equation with \m, b; t). The meaning of '~ 0' 
in (|62D is that the left-hand side of fl62j) which has the dimension of frequency must be much 
smaller than the frequency (energy) scale of the system. 

In view of (^T]) , we have the following expression for the adiabatic non-Abelian noncyclic 



geometric phase (|43|). 

= u£(t)r?(t) , (63) 
where w™ (t) is an TV x TV matrix with entries 

w n oab {t) := (n,a;R(0)\n,b;R(t)} . (64) 

If E n [R] is nondegenerate, TV = 1 and 

r:(t) = <(t)e^»W , where 7on (t) := / = / z(n; iZjdjw; i?> . (65) 



i?(0) J-R(O) 



The phase of f"(i), namely 



(n; R(0)\n; R(t)) 



1/2 



(66) 



(n; R(t)\n; R(0)) 

is precisely the Abelian adiabatic noncyclic geometric phase studied by de Polavieja and 
Sjoqvist M. 



8 Application: Spin 1 Quadrupole in a Precessing Mag- 
netic Field 

The simplest possible quantum system that allows for the occurrence of a non- Abelian geo- 
metric phase is a system with a three-dimensional Hilbert space and a dynamical invariant 
I(t) which has a nondegenerate and a degenerate eigenvalue |18| . If the system undergoes an 
adiabatic evolution, then the role of the invariant is essentially played by the Hamiltonian. 
In particular, the Hamiltonian must have a nondegenerate and a degenerate eigenvalue. The 
moduli space of all such Hamiltonians (for a three-dimensional Hilbert space) has the mani- 
fold structure of the projective space CP 2 , p3f. A thorough treatment of the problem of the 



adiabatic non- Abelian cyclic geometric phase for such a system is presented in Ref. In 
this section, we shall use the results of Ref. |28| to investigate the adiabatic non- Abelian non- 
cyclic geometric phase (t) for a spin 1 quadrupole interacting with a precessing magnetic or 
electric field. The problem of the non- Abelian geometric phase for a spin 3/2 quadrupole has 
been considered by Zee Mead |29| and Avron, Sadun, Segert, and Simon |3(J [31[| . For 
a bosonic system such as the spin 1 quadrupole considered here, one can show that Berry's 
connection one-form is a pure gauge and Berry's cyclic geometric phase angle vanishes . 
This result does not however generalize to the non- Abelian geometric phase, I jZSfl . 

Consider the quadrupole Stark Hamiltonian H = X(J ■ R) 2 , where A is a real coupling 
constant, J = (J\, J2, J3) is the angular momentum operator, and R = (R 1 , R 2 , R 3 ) is a 3- 
vector representing the magnetic (or the electric) field. For a spin 1 particle, this Hamiltonian 
has the form 

Xp 



H[R] 



( 1 + 2C 2 yfiCe-* 
V2 (e iip 2 



-2iip 



-vp 



\ 



-2iip 



V20 



-V2(e 
1 + 2C 2 J 



(67) 



where (p,<p,z) are the cyclindrical coordinates in the i?-space, i.e., 



p : = J(R}) 2 + (R?f , if := tan- 1 (R^R 1 ) , z := R 



( := z/p, and we have used the standard spin 1 representations of J,. 



In view of the general results of Ref. [28|, the eigenvalues of the Hamiltonian fl6T|) are 
given by 

Bk[R) = \p 2 ( 2 , and E 2 [R] = Ap 2 (l + ( 2 ) , (68) 

where R = (p,tp,Q. As seen from (0), if p ^ 0, then the Hamiltonian has two distinct 
eigenvalues. In this case, E\[R] is nondegenerate and E%[R] is doubly degenerate. A set of 
orthonormal eigenvectors of H [R] is given by [EH] 



\1;R) := iVf 1 



\2,1;R) := N^ 1 



( e-** \ 
v/2C 



V 



J 



( -y/2 Ce~ i,fi \ 



( 



V 



1 

o 



and \2,2;R) := (N^)- 1 



J 



-up 



\ 



V2C 
\ (1 + 2C 2 )e^ J 



(69) 



where N\ := y 2(1 + £ 2 ) and N 2 := y/1 + 2(' 2 . Note that these formulas are valid for p 7^ 0, 
i.e., everywhere except the i? 3 -axis. 

Again using the general results of Ref. or by direct calculation, we can compute 
Berry's connection one-forms A™. Doing the necessary algebra, we find 



Al[R]=d<p, and A 2 [R] 



pdip ve^dip 
ve~ tip dip adtp 



(70) 



where 



p :-- 



2C 



2c 2 



v 



a 



1 + 2C 2 1 + c 2 

C 



(l + 2C 2 )v / T+C 7 

1 + 2(1 + 2C 2 ) 2 

2(i+2c 2 )(i+c 2 : 



c(l - c 2 ) 
1 + c 2 ' 
l + 2(l + c 2 ) 2 
2(1 + c 2 ) 



(71) 
(72) 
(73) 



and c := C/(l + C 2 ) = z/Vp 2 + C 2 = R^/^R 1 ) 2 + (R 2 ) 2 + OR 3 ) 2 - Note that in the spherical 
coordinates (r,9,(p), we have 

c = cos 9 . (74) 



In view of Eqs. (|70|), Al[R] is a pure gauge. This was to be expected, for the system is 
bosonic [BTJ. The connection one-form v4„[-R] is not a pure gauge. In fact, unlike the spin 



3/2 systems considered in the literature |24|] even for the case of a precessing field where 9 
is kept fixed and ip varies, the adiabatic geometric phase associated with E 2 is non-Abelian. 

In the remainder of this section, we shall compute the adiabatic non-Abelian noncyclic 
geometric phase associated with the degenerate eigenvalue E 2 for a precessing field with 

9 = constant , tp = ipo + out , and <po, oj = constant . (75) 

First, let us consider the matrix 1^. We can write Eq. (|59"D as the matrix Schrodinger 
equation 

= %>)I» , (76) 

where 

K<p)d(p := -Al\R] . (77) 

Clearly, h((p) belongs to the Lie algebra of the unitary group U(2). In particular, it can be 
written in the form 

%>) = ~X>V (78) 
z 1=0 

where o"o stands for the unit 2x2 matrix, 01 with I e {1, 2, 3} are Pauli matrices, and 

r° := — {jj, + cr) , r 1 := — vcostp, 
r 2 :=vsm(p, and r 3 := a — [i . 

Substituting these equations in (^|), we obtain 

h(<f) = - [-(// + cr)(T - ^(cos (f o x - sin <p a 2 ) + {cr - n)a 3 ] , 

= ^/2[_( /1 + ff ) ff0 _ 1/ , 1 + (,_^] e ^/2 ) (79) 

where we have used the identity 



e 



-iipai/2 iipcTi/2 



3 



a ._ e vpa i /2 _ cog ^ ^ -f s in y3 V" eijk<Jk , for i ^ j . (80) 

fc=i 

In (|8H|), eijk stands for the totally antisymmetric Levi Civita symbol with €123 = 1. 

In view of Eq. (|79f ), h((p) is the Hamiltonian of a spin 1/2 magnetic dipole in a precessing 
magnetic field. Therefore, we can perform a unitary transformation of the Hilbert space 



to map it to a constant Hamiltonian. Under a ^-dependent unitary transformation 
U(<f), h({p) and T 2 (ip) transform according to |34j 



hUp) -> h'Up) = U((p)h(<p)U((p) ] - iU{p)^-U{p) ] 

dip 

Setting U(<p) = e"^ 3 / 2 in (P) and using (|79|), we find 



(81) 
(82) 



^' = 2 + (T ) (Jo - ^cri + (1 - A* + o")o"3] 



(83) 



For a precessing field, 9 and consequently c, /z, z/, and a are constant parameters. There- 
fore, h! is constant, and we have 



Substituting this equation in 



r^(v?) = e -ih'(r-v°) . 
we find 

a iip(J3 /2 — iti ( ip- tp ) „ -«V0 0-3 /2 



(84) 



(85) 



In view of Eqs. fl8"3] ) and (|55]), the matrix elements of r 2 (<yj) are given by 



r 2 

1 oil 



"/<-+J i< r - r -,,..-2 J cos [(^ _ (p )A/2] + i ^ ) sin[(^ - ^o) A/2] 



ol2 
2 

o21 



Z/ 



sin[(^ - <^o)A/2] 



(86) 
(87) 



r 2 

1 o22 



, , _^ e i(^)(^-vo)/2 e -i( v +^)/2 sin [^ _ ^ )A/2] , 
3 i(^ CT -i) (v -^)/2 / cos[( ^ _ ^ o)A/2] _ f ) Sin t(^ - V^o) A/2] 



where 



A := J(l + (J -ii) 2 + v 2 



l + 4c 2 (4 + 8c 2 + 7c 4 + c 6 ) 



59) 



(90) 



2(1 + c 2 ) 

Next we compute the entries of the matrix w 2 (t) of Eq. fl64l) . Using Eqs. Q69|) and doing 
the necessary algebra, we have 

/ 2c 2 















W ll2 


= w2 o21 










™o22 


= 1- 





1-e 



c(l - c 2 



= 1 _ ^(1 _ e -<(v-Vo)) ; 

1 - e-^-^) = - e - J (^)) , and 
' 2c 2 ^ 



1 + c 



1 + c 2 

[1 - cos(v9 - <p )] + i 



(91) 
(92) 



1 + c 2 



sin(y? - ip ) 



3 1 

1 + (<t + — + -)[!- cos(p - <po)} + in sin(v? - <p ) • 



(93) 



Having obtained r 2 and w%, we can use Eq. (|63| ) to calculate the non-Abelian noncyclic 
geometric phase f 2 = "U^r 2 - As seen from the above formulas the result will only depend 
on <p , if, and c = cos#. The expression for f 2 is rather lengthy Therefore, we shall instead 
give its trace IT 2 which is of physical importance, 



TT 2 = e *(M+o-)('^-¥'o)/2 [ % 



COS 



A((p - ip ) 



+ ysin 



A((p - ip ) 



where (94) 



X : = I e - i(< ^ o)/2 [6 + 7/i + 4a + (2 - 7/i - 4a) cos(y? - (p ) + 4z sin(<^ - (p )], and 

y := -L e -^-vo)/2( e i Vo _ ( 8j/ 2^ + e %+^)) + e i<p ^ 7/i _ 5) _ 3(2 + n) - 4(j 2 - 2] + 
8A ^ 

e^[/i(9/i - 19) + (14 - 13/i)a + 4a 2 + 10]} . 

One can check that for ip = <p , II 2 = 2, as expected. 

Furthermore, setting ip = y?o + 27r in (^j), we obtained the trace of the cyclic non-Abelian 
geometric phase (1591) which is given by 



n 



= _2e iff C^ +ff ) cos(ttA) (95) 

cyclic V 7 V 7 



In this case, we can easily compute the eigenvalues of the non-Abelian cyclic geometric phase. 
They turn out to be e ±in{A+1] . 

Now choosing <po — and c = cos 9 = 1/V3, which are the value used in Tycko's 



experiment fl3"2"], we have A = y / 889/24 « 1.24 and 



n 2 w ie ai0 ^{(2 + 4cos^ + 3isin^)cos(0.62^) + 

0.81[cos(^/2) + 2.42isin(^/2)]sin(v9/2)sin(0.62^)} . 

9 Conclusion 

In this article we used the theory of dynamical invariants of Lewis and Riesefeld to develop 
a general parameter space approach for the nonadiabatic geometric phase. We introduced 
a set of time-dependent gauge-invariant geometric quantities U n (t) which we identified with 
the trace of certain time-dependent gauge-covariant geometric quantities T n (t). For a cyclic 
evolution with period T, T n (T) yields the non-Abelian cyclic geometric phase. In the Abelian 
case, T n (t)/\T n (t)\ coincides with the Abelian noncyclic geometric phase studied in the lit- 
erature. Therefore, we identified the T n (t) as a non-Abelian noncyclic geometric phases 
factor. 



We also discussed the adiabatic limit f™(t) of t n (t). Again we showed that for the 
Abelian case r™(t)/\T™(t)\ is the known adiabatic noncyclic geometric phase. 

We have finally presented a through analysis of the Non-Abelian cyclic and noncyclic 
geometric phases for a spin 1 quadrupole in a precessing field. 

We wish to conclude this paper by pointing out the following observations. 

1. The original definition of the (nonadiabatic) non- Abelian cyclic geometric phase is 



due to Anandan ||20|| . As pointed out in Ref. Anandan's definition is identical 
with the one given in terms of the dynamical invariants. More specifically, the basis 
vectors \i(j a (t)) used by Anandan [p(J to yield the non- Abelian connection one-form are 
precisely the basis eigenvectors |A n , a; 9(t)) of our approach. In fact, our approach may 
be viewed as a means to identify Anandan's basis vectors \ip a (t)) with the eigenvectors 
of a Hermitian operator, namely a dynamical invariant I[9(t)]. 

2. In Refs. [J, Pati shows that the Abelian noncyclic geometric phase angle of Eq. ( pEBj ) 
may be written in the form 

In = I ^ , (96) 



c 

where fl n = i{Xn(t)\d\xn(t)) and \Xn(t)) is a properly scaled state vector. Although 
Pati terms Q n a connection one-form, he shows that indeed Q n is invariant under a 
gauge transformation. This is because Q is the difference of two connection one-forms, 
namely the connection one-form A n and another connection one-form P n which also 
depends on the initial state vector. As pointed out by one of the referees, it would be 
interesting to see whether Pati's results can be generalized to the non-Abelian case. 
Clearly, if one chooses a basis in the degeneracy subspace l H\ n {t) in which the non- 
Abelian noncyclic phase factor T"(t) is diagonal, then the diagonal elements may be 
treated as Abelian noncyclic phase factors and Pati's results may be used to express 
the corresponding phase angles in the form (|9~6"D provided that T n (t) (alternatively 
w n (t)) is invertable. In such a basis, f n = diag(gie^™, • • ■ , gi n e 1 ^™) =: G n e tSn , where 
'diag(---)' stands for a diagonal matrix with diagonal elements • • •, ge and 7^ are 
real, G n := diag(<7i, • • • , gi n ), and S n := diag(7^, • • • , 7^). In view of Pati's results, 
one may find an appropriate diagonal matrix of one-forms Cl n such that S n = fc& n - 



Furthermore, f2 n will have the form A n — P n for some matrix of one-forms P n . Now one 
may postulate that P n is a connection one-form, so that under a gauge transformation 
Cl n transforms covariantly. This together with the form of the diagonal elements of P n 
which is given by Pati || are sufficient to obtain Cl n in an arbitrary basis. A complete 
generalization of Pati's results to the non-Abelian case would require a set of properly 
scaled state vectors \x n , a; 9) satisfying Cl™ b = i(xn, a; 9\d\xm b] 9). The explicit form of 
the vectors \x n , 9(t)) is not known to the author. 

The function F introduced in section 7 may be used to relate the adiabatic and nona- 
diabatic Berry connection one-forms. Namely, the nonadiabatic connection one-form 
A n [9] is the pullback 36 1 of the adiabatic connection one-form provided 



that F is differentiable. This has been originally pointed out in Ref. ||19|| . However, in 
[[TIJ] the existence of F was assumed based on the evidence provided by the study of 
a magnetic dipole interacting with a precessing magnetic field. In the present article, 
we used the theory of dynamical invariants to establish the existence of F for a large 
class of quantum systems. In particular, it is not difficult to see that they exist for the 
systems possessing a dynamical group. This allows for the application of the holonomy 
interpretation of the nonadiabatic geometric phase using a fiber bundle which has the 
parameter space M. of the invariant as its base space [HJ. This bundle is the pull- 
back bundle -F*(A) of the U(Af) bundle A used in the holonomy interpretation of the 
adiabatic geometric phase [19|, |36|| . Note however that the function F depends on the 
adiabaticity parameter |27j . In general, there are certain values of the adiabaticity pa- 



rameter for which F becomes ill-defined or discontinuous. The above constructions are 
valid only for those values of the adiabaticity parameter for which F is a differentiable 
function. 

In our derivation of the noncyclic geometric phase, we assume that the Hamiltonian 
is a Hermitian operator. The generalization of our results to non-Hermitian Hamil- 
tonians is straightforward. One needs the machinery of the non-Hermitian dynamical 
invariants and their biorthonormal eigenvectors to obtain a non-Hermitian analogue of 
the noncyclic geometric phase. 
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